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ABSTRACT 


The objective of this project was to study the effect 
ciastaltedtrrtretal Cavrties on natural convection from a 
horizontal surface. Tests were run with Feron 113. Data 
for heat flux as a function of bulk temperature difference 
were carefully obtained. These data yielded Nusselt number 
as a function of Grashof number or Rayleigh number. All of 
these were then compared with the data obtained by Hiep 
(Ref. 2). Experimental results are presented for the heat 
meansrerstrom Horizontal circular disks, with and without 
meeiircial cavities. 

The artificial cavities were found to affect the 
mulral cCOnVection heat transfer from a horizontal surface: 
the heat transfer coefficients increased with the number of 


cavities. 
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Po ROU PON 


A. BACKGROUND 

iimeismoeen shown by Duncan [Ref. 1] that artificial 
cavities can have a significant effect on boiling heat 
meiiomenre if tne cavities are of the correct size it is 
possible to induce nucleation at lower superheats than would 
otherwise be necessary. Recognizing the enhanced heat trans- 
feorewiicn artificial cavities can induce in the boiling regime 
one is led to inquire as to the effect that these cavities 
might have on the natural convection heat transfer which 
precedes boiling. Hiep [Ref. 2] undertook an experimental 
mivestigation of the effect of artificial cavities on the 
fecural convection heat transfer from a horizontal circular 
disc. Tests were run with water and Freon 1173. For water 
Mmemonnior ket. 3] ran tests with a mirror finished test 
Surface without the artificial cavities and plotted the data 
in terms of dimensionless parameters, Nusselt number versus 
Grashof number. Hiep started with seven artificial mercies 
for water, and gave the correlation result of Nusselt number 


in term of Rayleigh number as: 


Nu = 0.38 (Ra)?°*? (1) 
For Freon 113 Hiep ran the tests with a blank and seven- 
hole test sections, again he gave the results in terms of 


Nusselt number versus Grashof number and Nusselt number against 


im 





Rayleigh number the correlation results for blank test section 


Was 
iy 2 0.48 Upro oes (2) 


The experimental apparatus that O'Connor [Ref. 3] and 
Hiep [Ref. 2] used was almost the same. Test surface had a 
diameter of 0.895 inches, installed in the center of a 
cylindrical tank 9.00 inches in diameter and 4.5 inches high 
and set at the same level as the tank floor (see Figure 3). 

A controllable resistance heater as a heat source was housed 
in a copper block attached to the bottom of test section. 
mecaliiprated resistor in series with the heater was used to 
determine accurately the input power to the heater. Eight 
thermocouples, wired in series were used to measure the 
average bulk temperature of the fluid and one removable 
thermocouple was used to provide a check on this temperature. 
To record the temperature distribution in the stainless steel 
cylinder, four stainless steel sheathed thermocouples were 
equally spaced along the axis of the test section. The 
temperatures at these locations were obtained using a Dana 
digital voltmeter. 

Since all the results and discussion of Hiep's work 
depend on the accuracy of determination of the test surface 
and fluid temperature the calibration of all thermocouples is 
extremely important. Hiep used the standard thermocouple 


tables instead of the average of his calibration value. The 
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uncertainty value in Hiep's thermocouple reading was 0.5°F, 
waemeresults in a large uncertainty for heat transfer 
coefficient. The change in Nusselt number between the blank 
test section and the seven hole test section as determined by 
moo icewithin the experimental uncertainty assigned to the 
Nusselt number. 

The uncertainty in Hiep's experimental data, therefore, 
meds tO make it difficult to draw a firm conclusion as to 
the effect which artificial cavities have on the natural 
Gomvection heat transfer process. His experiments were 
limited to two test surfaces, blank and seven cavities. He 
also jlimited the range of surface temperature to those below 
the saturation temperature of Freon 113. In other words the 
experiment did not go into the boiling phase that Duncan 


[Ref. 1] obtained. 


Be filesolS OBJECTIVES 
The previous discussion has led to the objection of this 
experimental works, that is: 
ie To investigate the effects of the drilled cavities 
Cmematiiral Gonvection heat transfer in Freon 113 
with greater number of test surfaces than Hiep 
obtained. 


ae fjemontain results for the boiling phase of Freon 113 
to see if it agrees with the Duncan results. 


oe Try to keep the uncertainty value of heat transfer 
coefficient and other dimensionless parameters as 


low as possible. 


LR: 





TI. EXPERIMENTAL APPARATUS 


A. DESIGN CONSIDERATION 

ile The experimental apparatus used in this study is 
the same as had been used by Hiep [Ref. 1]. A complete 
@esSeription can be found in his thesis and the thesis of 
O'Conner [Ref. 3]. A schematic diagram of the experimental 
set-up is shown in Figure 3. For convenience the ice-bath 
used by Hiep was replaced by an automatic electronic ice-point 
for the reference thermocouple. This automatic ice junction 
was set to 32°F as reference temperature. A Dana digital 
voltmeter that had been used to read and record the thermo- 
couples was replaced by a Hewlett Packard digital voltmeter. 

Za AeFreon condenser was also added to recover the 
Freon that evaporates during the experiment. This condenser 
was connected to the top plate of the tank as shown in Figure 


eemond was cooled by tap water. 


B. fits oECTION-HEATER ASSEMBLY 

Three test sections were used in this experimental work: 
one blank, one with seven holes which had been used by Hiep 
and one with thirteen holes (cylindrical drilled cavities). 
All three test surfaces had the same diameter, that is .896 
inch. Four 0.040 inch diameter holes were drilled along the 
Begin O17 each test piece to a depth of 0.375 inches, separated 
axially by about 0.5 inch and 90° apart from each other. 


These holes were for sheathed thermocouples installation in 
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Figure 1. 
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Figure 5. Photograph of Thirteen-Hole Test 
Section After the Test. 
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each test section as shown in Figure 6-8. The thirteen-hole 
test section had thirteen cylindrical holes of the same size 
as the holes drilled on seven-hole test section, that is 
0.0375 inch in depth and 0.015 inch in diameter, drilled as 
shown in Figure 8. The heater section [Fig. 9] was the same 


as used by Hiep [Ref. 2]. 


e. INSTRUMENTATION 
2 Temperature Measurement 

ISA type T copper-constantan thermocouples were used 
to measure the temperature distribution in stainless steel 
meee section in order to get as accurate as possible in 
temperature distribution. These four thermocouples were 
Siemacned in Stainless steel! (sheathed diameter = 0.035 inch) 
and the temperature of the surface was found by using a linear 
extrapolating method (see Appendix II, Fig. 18). 

The bulk fluid temperature was determined by using 
eight copper-constantan thermocouples wired in series (four 
pairs) manufactured using a Dynatech welder. These thermo- 
couples were wired in series into an insulated switching box 
so the reading gave an average temperature for four pairs of 
these thermocouples. A single removable thermocouple was 
used to provide a check on this temperature. The sheathed 
thermocouple leads were securely fastened to the bottom 
phenolic plate and then to a Jones bar arrangement. Because 
the sheathed thermocouples were difficult to manipulate and 


Subject to twisting, and thus shorting, giving erroneous 
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reading, the more flexible wires were used for the bulk 
temperature thermocouples were used as extension wire into 
the switching box. A reference thermocouple was inserted 
-jnto the automatic ice-point reference. A Hewlett Packard 
Gigitai voltmeter accurate to 0.001 millivolt was used to 
read the output of the thermocouple in millivoit. 

Zi Power to Heater 

A Lambda regulated power supply was used to provide 

input voltage to the heater. A calibrated resistor was placed 
in series with the power supply and heater. This resistor 
was mounted on a large heat sink and was isothermal and 
tnerefore of constant resistance for all of the experimental 
runs. With this arrangement, a voltage reading across the 
resistor and across the heater were recorded. Knowing the 
fomeeage, and resistance, the current through the circuit could 


be determined by Kirchoff current law: 


Ty = Ip = Vp, less) 
where I, = the current across the heater. 
Ip —ithemelmmenmt —decross Une resistor. 


R = 2.013 ohn. 


Vp = voltage across the resistor. 


and the product of this current and voltage across the heater 
gave the input power, i.e., 
ee 


or 





mente oie could be aecurately reset for different runs. 
The voltage across the heater and across the calibrated 
resistor were determined using the Hewlett-Packard digital 
voltmeter, with an accuracy of 0.001 volts. This power was 
not used to determined the heat flux but only to set the 


Smerating point from run to run. 
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i eemmnenNAL PROCEDURE 


AY. fees?) SURFACE PREPARATION 

It was very important that the test surface be free from 
any machine marks and scratches that could influence the heat 
transfer other than the artificial holes. For this reason 
the surface was mirror finished. The machine marks were 
removed by using a Buehler Metallurgical Manual grinder; both 
O emery polishing paper and 3/0 emery polishing paper were 
used. After that, the surface was washed with detergent and 
warm water to remove the previous abrasive and then washed 
with alcohol. In order to get a mirror finish the test 
Surface was then wet polished on a metallurgical wheel covered 
with velvet. The wheel was impregnated with one micron diamond 
dust. The test surface was placed face down on the wheel and 
a low speed of rotation was used. During this work alcohol 
was used to wet the surface at all times. After completion 
of polishing, the test surface was washed with alcohol again 
and was dried with a dryer, and then wrapped with special 
paper to prevent contamination and touching. Before installa- 
tion the seven-hole and thirteen-hole test sections were put 
in an ultrasonic cleaner for ten minutes, were washed with 
distilled water, rinsed with acetone and dried. 

The seven-hole and thirteen-hole test surfaces were 
drilled before the mirror finish, by using a Sphinx spirec 


pivot drill, 0.015 inch diameter. The cavities depth were 


2] 





controlled by using a feeler gauge on the lathe. A ratio of 
approximately 2.5 was desired for cavity depth to diameter, 
SO the depth was about 0.040 inch; the location of seven 
holes and thirteen holes were as shown in Figures 7 and 8 
respectively. 

A height gauge which had a scale readable up to 0.0001 
inch was used to determine the location of the center of the 
thermocouple holes relative to the test surface. This gauge 


meeemarsoO Used to measure the diameter of the test surface. 


B. APPARATUS ASSEMBLY 

The main parts of the apparatus were the test section, 
heater assembly, and the tank assembly. The test section- 
heater assembly was securely bolted to the phenolic base 
plate, as described by Heip [Ref. 2] and once the insulation 
was in place and the can installed, the phenolic plate that 
Served as the bottom of the tank was carefully fitted to the 
test section. The level of the test section was adjusted 
from the bottom base of the assembly by screws. The phenolic 
bottom plate of the tank was secured by three threaded rods 
120° apart, with the phenolic base plate. These three rods 
were also used to adjust the test surface in level with the 
premolic bottom plate. For the tank assembly four "0" rings 
of 9.0 inch diameter were used for sealing to make sure that 
the cylindrical enclosure tank was leak proof. Two "0" rings 


for each side (top and bottom) of the tank were used. The 
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eemarical plexiglass container was then put in place and 
the top installed, six threaded rods, 60° apart were used to 
secure the top plate to the bottom plate of the tank. The 
glass laboratory condenser was installed on top of the test 
tank to recover any vaporized Freon. Two of the threaded 
rods were 18 inches long so they could be used to hold the 
condenser and insure that it would sit tightly on the 1 inch 
hole at the top plate of the tank. This 1 inch hole was used 
memour the fluid into or out of the tank. To hold the 
condenser an adjustable holder was connected from each of 
the 18 inch rods to make the condenser stand still. The 
condenser was held tightly on the top but was also removable 
to allow filling of the tank. 

Six small holes were drilled through the plexiglass top 
of the tank to insert four pairs of thermocouples in the 
tank to measure the average bulk temperature of the fluid, 
for each pair one of the thermocouples was put near the 
bottom of the tank and the other one was placed near the top 
so that the average of fluid temperature all over the tank 
could be measured. Two holes, one in the middle of top plate 
aigeecne Other on the side of top plate, were used to insert 
one removable thermocouple to the middle of the tank to 
provide a check of the bulk temperature read by the other 


four pairs. 


29 





€. TEST PROCEDURE 
AP enter aburne 

A complete data run consisted of starting at the 
lowest heater setting, waiting until steady state was reached 
(usually about ten hours for each data point). The data was 
recorded every hour during the daytime, and at longer inter- 
vals during the night. After steady state was reached and 
the data for that power setting was recorded, the power to 
the heater was increased by approximately 5 watts. This 
procedure was repeated until boiling occurred from the surface. 
The power was then decreased in 5 watt increments and the data 
was recorded after steady state was reached. This procedure 
was repeated until the starting Bownt was reached. | The 
procedure of taking data during both a heating and cooling 
process was done to check on any hysteresis effects which 
eont exist. 

The data recorded for each setting were the bulk 
temperature, the temperature of the individual sheathed 
thermocouples in the test section, and the voltages across 
heater and calibrated resistor. When the ene temperature 
and the temperature of the test section changed by less than 
half degree F per hour, steady state was assumed to be 


reached. 


s 





lamp Determination of Nusselt, Grashof and Rayleigh 


Numbers 

For each run, at al] input settings, the temperature 
distribution in the stainless stee] test section are plotted 
and found to be linear. The surface temperature could be 
determined within 0.5°F from the plot of temperature versus 
thermocouple location by extrapolation, as shown in Appendix 
B. The properties of the fluid were determined at the film 
temperature which is the average of the bulk temperature 


a) and the test surface temperature (To) tee 


The thermal conductivity of 304 stainless steel 
(Keo) was calculated using the data from Touloukian [Ref. 10] 
at the average of the temperature of the two nearest thermo- 


couples to the test surface i.e., 


This corresponds to the average temperature of the 
length 6x. A simple energy balance at the surface using 


Fourier's Law of Cooling gives: 


61. 
meee 2 Dict.) = het, 
K aor) 
oS 2 O 
h z “Ss (GTS) Buy Die — tit =F (5) 


3] 





———<—-* <<” —* 


The Nusselt number was obtained from the definition: 


Nu = t. (6) 


psc witliamamevalwes ot ti trom (3) into (4) gives: 


Nu = ss 2 | mn OP) 
] i 
The Grashof number was calculated using the diameter 


mare test surface as the characteristic length: 


ae SD (8) 


The Rayleigh number is the product of the Grashof 


number and Prand!it number: 


Ra = Gr-Pr (9) 
where 
Py “ a . ylp Caos 
—_,: Ke 
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Re Vi eee oCUSS FON 


It was found that a minimum of ten hours was necessary 
for any input setting to reach steady state. Once steady 
state was reached, very good agreement between Nusselt number, 
Grashof number and Rayleigh number, based on test surface 
diameter, obtained during the heating and cooling phase of 
the run resulted. The uncertainty bands are shown on the plot 
Sueeeach curve. The investigation for the blank test section 


=O to 350x10°. for 


and 25x10° ee 


covers a range of Rayleigh number from 30x10 


6 6 


seven-hole test section from 50x10° to 250x110 


350x10° 


for thirteen-hole test section. A Summary of the 
meemrits for various data runs are shown in Tables I, II, and 
Meee FOr each data run the Prandtl number, Nusselt number, 
Grashof number and Rayleigh number were calculated. These 
dimensionless parameters were determined using the test 
section surface diameter as the characteristic length. As 
described in the experimental procedure, data were taken at 
each setting as the apparatus was heated with increasing input 
and then again as the inputs were decreased, cooling the test 
surface. The cooling phase data are denoted in the table. 
Surface temperature of the test section was determined 
by linear extrapolation from all four sheathed thermocouple 


reading. This gave the surface temperature with an error of 


approximately .5°F.° 


oS 





"aseud ButL{ood ay. paeqzussaudauy, :d330N 
8L°OE€ 655°9 O29°E 976°L 6S€°88 bao €°S8 mje || |h 
6S°EL €S6°9 vL6°8 SZ 9/8° LLL [oe 6 o6 ¥VEZE 
G9° POL Gila Z le 2 Lt [30 Oi Oral ilk Gaull x6 0° P 

0°022 pov ll eA [ok Z 906°VtL 9°€8 ae hall x2L°6 
Eo Se 90¢°Z SO Gis ScG@.e Sel eS €° 98 Se lt ~9GE°LL 
POLE Gilgiae ee GLO°?Ss Git ieiee Gomeoed 0°60L c'?7Sl GEL ve 
Oise ie) Ti p9° LP SS Ogee Ly SiO 0°6tl GS9°62 
iar G29'°L AS PA 88S °¢ oy Aan 23 éc0il 8° Trl OS Pe 
oe ale Coy 1 GSO'LE pie Coac Cae | 0°88 ve 2c 5137 (LIL o 
Cue LO Cvyy PSG 88L°¢ Sea ait L6°98 Cmaa. 9€5°6 
CL° SSL 68€°Z DS) a ONS) [LZ LS“ Lvl 6°08 9°SLl py°8 
5) SVE Ih O9C°L LO 1c Oileac VO. ool 9°€8 8° PLl Ov9'L 
ye 10h 1G) ae Vowel v80°¢ V6) ileal 0°L8 Sc" OoL COGeL 
91°29 916 °9 6€2°8 Gu0l-< OO el Oe Als Opis 
Ob €e£g°9 Cree Cee | GZ1°06 Gar bs Loin 
g-OLXRY (99) g-OLX45 (ON) nN q, $ Shout 


udL19aS 32Sa] yuel,g vos SatNsay so Auegwwns 


Dassieey. 








TES. 
6¢ 601 
Beet 
8°80¢ 
SSI 
nec 
Ov v9l 
09° 801 
eS) 


g-OLXEU 


SiGe 


> 
bO7 


UOL}ZISS YS9] S9LOH-UdA|DS UOS SZ2[NSAY JO Auewwns 


Ore 9 
Uae el 
OUT: id 
026 Le 
GS) ss 
Gee Oh 
Os Ieee 
Uc See 
eo © 


9 


BOK) 


vv0'¢ 
Sie ec 
eae ae 
GC a? 
5) SP eee 
cO¢c ¢ 
LO leg 


OiCie ee 


0) a2 


(ON) 
B07 


el alelaiy alt 


"aseyud Bultpoods ayy 


Space 
Lio She 7: 
cL 6tl 
L0°p9l 
OE eG 
S102 GgiL 
ole SIF 
spebel G'S |b 
08° SOL 


My 


0°64 
GL 
Cs 
bp L8 
L° L6 
0°98 
C78 
Viseow 


PaqyUuasoudau , 


0°06 
3) el) (L 
O°PLl 
Orrell 
Spey 08 IL 
IL IL 
o eu 
8° O0L 
L°88 


Sy 


>320N 


SCV ¢ 
x66¢ § 
*96L L 
*LbL Ol 
Oo vA 
Sey OH 
Gree 
Ones 
Cay ec 


(322M) 
ynduqJ 


oe 





Crees 
UZ ce 
Pemelt 
09°88 
09 Ord 
OL°8¢ 
OL me0e 
OD 062 
LiL 16) 
Ol ol 
He col 
20 eal 
OS Pe 


9 


_OLXed 


OES) 33, 
876° 9 
Sore 
Sie. 2h 
og 
Gore 
1 Sie) % 7 
GS) Ie 
EG 2 
Shsitea ae! 
Wis) [La 
Seley 8) 
We 8) 


(49 ) 
Boy 


UOLI99S 3S9] SLOH-usaqulLuL wos Sdinsay jo Auewwns 


S07 
018 °8 
Oe al 
G66 °?¢ 
LpD Eee 
U6 00 
O8l° SP 
OD ey 
Nene) as 
O92 Te 
Cet 
Elka 
Gee 


9 


_OLXAD 


a NABI 
Cicada ec 
Gail 
SVG 
D8 ¢ 
CVG 2 
SUC 
Oe 
A 
WEG 
O8l°¢ 
Coelene 
ie eiloha aL 


(ON) 
boq 


Lava 


"aseud BuL,ood ay. payUassudau xy 


60°OLL 
Ob SEL 
99 ee V Sul 
V Le SOCn 
Dee Oe 
OS) 0 ca 
eySAOl irgs) 

08° 99D 
Ops’ Lle 
Soe ios) 
Ose SE 
(ASS S) © evra 
Cuca IG 

ON 


9°08 
L°08 
c° 18 
Lac 
L S16 
GS OOL 
ean 
Perel 
Coen 
Cuscic 
Goo 
oy A 


€° 88 
b' 6 
St 10h 
G 20ne 
aieCel 
Scull 
i (Sih 
oat 
GeO cil 
Ck 
LC 20k 
86.6 
p28 


aa 


:O2ON 


¥8¢2G°L 


¥PE°€ 
¥8L°S 
x62 6 
xl? 9l 
x60 9¢ 
ao We 
Sida g¢ 
er gal 
oy a an) 
SS LS 
MAES te 
eval 


(322M) 
4ndu] 


36 





A log-log plot of Nusselt number versus Grashof number 
and Rayleigh number based on test surface diameter is shown 
momeeacn test Surface in Figures 10 to 17. The curves show 
that the results fall into two regimes: 

lis: Natural convection which is the purpose of this 
investigation. The Rayleigh number for this range falls in 


6 to 180x10° FOr venese bnree test “striace results. 


between 30x10 
The curves show the straight line and almost have the same 
slope (see Figures 16 and 17). The correlation which results 
for the natural convection range of each test surface with 


Freon 113 as the fluid are: 


Nu = 0.48(Ra)9*29 30x10°<ra<180x10° for blank test 
Secrron: 
Ocoee = a0 1022 ka180x10° for seven-hole 
test section. 
Nu = WALD kanes 20 OMERe <1 ose for thirteen- 
hole test section. 
ae Boiling phase, the results of which must be 


considered with some caution. This is because the bulk 
temperature of the fluid was always below saturation and 
therefore the boiling washighly sub-cooled. In addition when 
mucleation occurred it began at the edge of the test surface 
in the space between the test section and the phenolic plate. 
Subsequently nucleation did take place on the test surface 

as can be seen by the cavitation damage visible in Figures 

4 and 5. The Rayleigh number for this range started from 


200x10° and up. The curve in this boiling phase is shown as 


3/7 





a_straight line with slope greater than the natural convection 
line (also see Figure 16 and 17). 

Miele natural convection range the results show that the 
meametranster coefficient of the fluid increases. This is 
consistant with Duncan's result but disagrees with Hiep's. 
Serpe stated that as the number of drilled cavities increase 
meme neat transfer coefficient decrease. The discrepancy with 
meeps result can be attributed to several factors: 

ie. The uncertainty in thermocouple calibration which 
mep used was 0.5°F which led to a large uncertainty value in 
Surrace temperature. Surface temperature was determined by 
linear extrapolation from two sheathed thermocouples near the 
Surface. This has a large effect on the results which depend 
on the surface temperature. For the present experiment the 
value of the uncertainty in thermocouple calibration was 
reduced by carefully calibrating the thermocouples and uSing 
the average value of the calibration. This reduced the 
uncertainty in the thermocouple measurement to 0.2°F, which 
results in a smaller error in heat transfer coefficient. 

a. The number of test surfaces that Hiep used (blank 
and seven-hole) was insufficient to conclude that heat 
transfer coefficient decreases while the number of cavities 
ieereases. that is Hiep’s results for the seven-hole test 
section fell within the uncertainty band of the blank test 
section. The present experimental work used a thirteen-hole 
meee SECTION in addition to the seven-hole test section. 
These results show a consistent trend of increasing Nusselt 


number with increasing number of cavities. 
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V. CONCLUSION 


The results of this experimental work lead to the 
mom lOWwIng conclusions: 

i The artificial cylindrical drilled cavities affect 
the natural convection by increasing the heat transfer 
coefficient of fluid as shown in Figure 16 for the log-log 
plot of Nusselt number versus Grashof number and Figure 17 
for the log-log plot of Nusselt number versus Rayleigh number. 

Ze The correlation results for this experimental work 


mene horizontal circular heated disk in cylindrical enclosure) 


are: 
Nu = O.48(Ra) °°. TOr blank tect esect ion 
Nu = VSatme coe, for seven-hole test section. 
Nu = Bae for thirteen-hole test section. 


All of these correlations were in the range of Rayleigh 


6 6 


number from 30x10” to 180x10°. For the value of Rayleigh 


6 


number from 200x10° to 350x10° the slope 1S increased as seen 


in Figures 10-17. 
Se. The log-log plot of Nusselt number versus Grashof 
number and Nusselt number versus Rayleigh number shows that 


the curves were straight lines and had nearly the same slope 


6 6 


in the range of Grashof number from 3xl0° to 20x19” and 


Rayleigh number from 30x10° to 180x10°, Between the range of 


6 6 


Rayleigh number from 180x10° to 200x10° the slope started to 


47 





change for each curve, and then the curves were straight lines 
again for Rayleigh number greater than 200x109 as shown in 
Praures 16 and 1/. 

4. Determining the surface temperature by the method 
of linear extrapolation as shown in Figure 18, gives the 
meeres Within an error of 0.5°F. This is insufficient 
because the results show that if the surface temperature is 
emanged by only .2-.3°F, the value of Nusselt number varies 
by about + 10% in the lower power input setting. This made 


the log of Nusselt number change by about +122. 
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APPENDIX A 


THERMOCOUPLE CALIBRATION PROCEDURE 


The accurate determination of the test section and fluid 
temperatures waS a major part of this experiment since all the 
results and discussion are based on these temperatures. For 
this reason, precise calibration of the sheathed and bulk 
temperature thermocouples was necessary. When the calibration 
was accomplished all the wiring and recording instruments that 
would be used in the experiment’ were not changed. 

A Rosemount Calibration System, with a constant tempera- 
ture oil bath, was used for the calibration. The ice junction 
thermocouple which was sealed in a glass tube filled with oi] 
was placed in the electronic ice point reference junction. 

The four sheathed thermocouples and nine bulk temperature 
thermocouples were suspended several inches into the oil bath. 
A Platinum Resistance Thermometer in conjunction with a 
commutation bridge was used as a standard. The calibration 
was conducted over a range from 80°F to 220°F and a maximum 
error for the thermometer for this range of temperature was 
0.004 millivolt. 

The results of these thermocouples reading compared with 
standard thermocouple tables. The deviation from the table 
values ranged from +.004 to +.016 millivolt for sheathed 
thermocouple, and the deviation from table value ranged from 


+.006 to +.016 millivolt for bulk temperature thermocouples. 
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Since all the deviation is in the position range the decision 
was made to use the reading value and subtract the average 
difference from the table value. The values of uncertainty 
for sheathed thermocouples was 0.2°F and uncertainty of 0.2°F 


for the bulk temperature. 
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APPENDIX B 
SAMPLE CALCULATIONS 


a 


For Thirteen-Hole Test Section: 


Input Power = 3.382 watts 
(Data recorded 6 October 1975) 


A. ee RMINATION OF THE SURFACE TEMPERATURE 

A graphical method was used to determine the surface 
Mmemoerature. For each input setting a plot of the temperature 
versus the thermocouple distance from the surface, as shown 
in Figure 18, was made using all four reading of sheathed 
thermocouples. The plots were linear and the surface tempera- 


mee Could be obtained within an error of 0.5°F. 


B. CALCULATION OF THE PRANDLE NUMBER 


From the definition: 


C 6 
pr= Y=. pb = yh (10) 


where fluid properties were evaluated at film temperature. 


ee | 


[T, = aa 5 bD = 86.6°F] 


Fiuid properties at Ty, = 86.6°F [data from Table I] from 


Thermophysical Properties of Refrigerants [Ref. 4]]. 


24 





Power Input 3.382 watts 
Bulk Temperature 79,4°F 


eis 


100 





90 
0 ae a4 (6 7S T3021 de eG TC ee) 
Distance from Test Surface, Inches. 


Figure 18. Determination of Surface Temperature. 
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y= «61.519 1b, /ft-hr 

ke = .04277 BTU/hr-ft-°F 
o- _° 

, 230 BTU/1b_-°F 


Using equation (10), we have 


pee clei) (9 2300)) 


c. CALCULATION OF THE NUSSELT NUMBER 


Using the equations: 


eeenGili= ) 
-— SS ie (5) 
(6x) (8T,) 
and 
io = (6) 


th 


The value of Keg was determined for an average tempera- 
ture of cylinder over the length 6x. The average of the two 
nearest temperatures to the test surface in the stainless 
Steel cylinder were used. The value of Keo was obtained from 
a plot of temperature versus Kec based on values from 
Touloukian [Ref. 5]. Again fluid properties were evaluated 
at the film temperature using the data from Table II and 
Table 13 from Thermophysical Properties of Refrigerant 113 


[Ref. 4] and the following data were obtained for 3.382 watts. 
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T, = Ooo 


ime = 100.9% 

Keo = 8.926 BTU/hr-ft-°F 
ee eG aa 

6T, = 14.6°F 

sT, = 5.00°F 

6X - CU ye Saete ns 

; = 96.6 Ib /ft.9 


The value of 8 were calculated by definition from the 


thermo-dynamic properties [Ref. 6] 





1 do eel 0 
B= =o = = — {-_-) 
Cte 0 Ty To : 
where 
P1*Pg 
as 
T,-Tp) = Temperature difference of the given value 


p> and PQ: RmOMm Ici ey I oamincema. cae 


For this run, we have 


4] 


B = 8.668x10— oF 


Using equation (5) and (6), the Nusselt number is obtained: 


Beg26( 5.) 


: Zo 
h = “HaT7S 14.4) 74.2395 Bane tee 


_ 74.235(0.07467) _ 
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De CALCULATION OF THE GRASHOF NUMBER 


Using equation: 


3 st 


ee (g8/ 2) D 


The Grashof number is obtained 


oo a SEES: BGh EO ne 
YY = en 


(1.519 )* 
56.9x3600 


(0.07467)°(14.4) 


6 


Ce ex 10 


an PyecULATITON OF THE RAYLEIGH NUMBER 


Using equation: 


Ra = Gr-Pr 
aus 
ae = (8 CN) ea 
= 7.202x10/ 


He, 





APPENDIX C 


UNCERTAINTY ANALYSIS 


In this experimental work the method proposed by Kline 
mmeemicClintock [Ref. 7] was used to obtain the uncertainties. 
Several assumptions and approximations had to be made in 
order to develop a reasonable estimate of the uncertainties 
for each variable. The following table is summary of these 


quantities. 


TABLE IV 


Uncertainty of Variables 


Tales (Ref. 4] 
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Variable Basis for Value Uncertainty (U) 

Cp assumption based on table Src. y SW Vee 
values [Ref. 4] 

D accuracy of measurement Beit ite 

Ke assumption based on table 1x107° BTU/hr-ft-°F 
values [Ref. 4] 

Keg assumption based on table 5% 
Values [PRef. 5] 

To based on Scale of GGannl calm o ar 
method 

Te thermocouple calibration Oe 

Ty 2712 thermocouple calibration Oise or 

6X accuracy of measurement 4.2x107" fete: 

B assumption based on table Ise yk 
Values [Ref. 4] 

: assumption based on table _1x10°* 1b /#t.° 
values [Ref. 4] 

u assumption based on table 2x10? Ne eis 





The method that proposed by Kline and McClintock 
[Ref. 7] for determining the uncertainty values is the 
Second-power equation. 

Uncertainties in the temperature differences 61, and 


615 were calculated from 


2 2% 
U = U SU Us 
a eS Vis Us 
ale 
: 2 2,° ° 
Sl Oe Sy Oe) re oilers 
and 
Ver erst, = E(Up ye + (Uy 21" 
2 e3 2 e 
2 2? 
SO 2 i) ai Oi?) ee) eer 
Where Us = uncertainty in T. 
S 
sil = uncertainty in Vie 
Same as U U = uncertainty in T, and T, respectively. 


Uncertainty in the convective heat transfer coefficient 


(h) was determined by using 


Maevalues of uncertainties are obtained from Table IV. 


S17) 





The Nusselt number jis defined as 


= 
CO 


Nu = 


Giz) 


~~ 
—h 


and used of the second power equation for uncertainty in this 


number becomes 


2 2 vu, 27 
Muy 7 EC) # GD + (et I (13) 
tee Uy 
D ii 
Tt can be shown that the values of (>) and qe re 
fF 


very smal] and can be neglected. Therefore, 


U 


Nu _ 
Nu" h (14) 


To find the uncertainty of the Grashof number use this 





relation: 
Gr = (gare) p? 6T, 
SO 
Mar 2 (oy + (Up) + (24) + (34D) nyty (15) 
r B p 61, 


The values of U are found from Table IV. 


And by definition for Prandtl number is: 


HG 
Doo pl ae 


K ¢ 
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So 


2 2 U 2 4% 


U U U k 
pee LE) + (eR) + (GDI (16) 


i 
against all the values of U are found from Table IV. Find 


the uncertainty value for the Rayleigh number, 


where 
Ra = Gr-Pr 
2 2 4 
U U U : 
eet CE) + RA (17) 


A summary results of this analysis are given for each 
mene in fable V for blank test section and in Table VI for 


meven-hole test section and for thirteen-hole in Table VII. 


og 





Power In 
(Watts 


- 


Uncertainties in the Experiment (Blank) 


Nusselt 
Number 


AO 
oe 
Ve 


09 
[2 
Ce. 


se 
70) 
nO 
As): 
. 30 
28 
eo 


TABLE V 


Percent Uncertainty 


60 


Grashof 
Number 


8.40 
4.04 
Sek) 
Ae) 
“NG 
@ (0s 
Be che, 
sss 
or 
Eu 


Rayleigh 


Number 





TABLE VI 


Uncertainties in the Experiment (Seven-Hole) 


Power Input Percent Uncertainty 
(Watts) Nusselt Grashof Rayleigh 

Number Number Number 
2.429 hast 606 5.64 
5. 300 Ces 0) 2.90 Sa74 
me / 9) 5.98 oro ou oO 
ime. 060 5 016 2.09 3.19 
14.970 oye 1.95 Silke 


6] 





TABLE VII 


Uncertainties in the Experiment (Thirteen-Hole} 


Power Input Percent Uncertainty 
(Watts) Nusselt Grashof Rayleigh 
Number Number Number 
e167 23.20 10.26 Ils 
a a opey 8.48 4.00 4.64 
pe 758 6.42 (age 3.64 
9.176 5.66 (ae |S Seed 
6. 300 oo i936 3.14 
26.280 oon i 210 Silks 
ea. 85 S48 94 Shin ING 
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